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Abstract 

In [3] the remaining area of the left-lower convex hull of a Poisson point process with intensity 
one in the first quadrant of the plane was analyzed, using the methods of [5] , giving formulas for 
the expectation and variance of the remaining area for a finite interval of slopes of the boundary 
of the convex hull. However, the time inversion argument of [5] was not correctly applied in [3], 
leading to an incorrect scaling constant for the variance. The purpose of this note is to show 
how the correct application of the time inversion argument gives the right expression, which is 
in accordance with results in [9] and [1]. 

1 Introduction 

Let P be a homogeneous Poisson point process of intensity 1 in the first quadrant of the plane. We 
consider the left-lower convex hull of this point process, as in [5], where a Markovian jump process 
of points W{a) is introduced, by which one can "walk" through the vertices of the convex hull. The 
points W{a) are defined in the following way. 

Definition 1 For each a > 0, W{a) = {U{a),V{a)) is the point of the realization of the Poisson 
process V on such that all points of the realization lie to the right of the line of the line x+ay = c 
which passes through W{a). If there are several of such points (which happens with probability 
zero for fixed a), we take the point with the smallest y-coordinate. 

In [7] the following result is proved for the joint distribution of the number of jumps and the 
remaining area, corresponding to a (slopes) interval [a,b]. 

Theorem 1 Let N{a,b) be the number of jumps in the interval [a,b] of the process W, as defined 
in Definition 1, and let D(a, b) be the area of the union of the triangles Ti, corresponding to points 
of jump Oi £ [a,b], as defined in Theorem 2.1 of [7]. Then: 



(A log(6/a)) (iV(a, b) - i log(5/a), Z)(a, b) - i log(6/a)) A iV(0, S), b/a 



oo, 



where N{0, S) is a bivariate normal distribution with expectation and covariance matrix defined 
by 
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This result is used in [7] to prove the unpublished result in [9], which gives a central limit 
theorem for the joint distribution of the number of vertices and the remaining area of the convex 
hull of a uniform sample of points from the interior of a convex polygon, see Theorem 1.1 in [7]. 
This extends the results in [3] and [5] . 

Theorem 1 is illustrated in [7] by a simulation of the jump process, generating the vertices and 
the areas, without first having to generate the points of the Poisson process for which one considers 
the convex hull. In this simulation Theorem 2 below is used. It is the purpose of this note to 
provide a proof of this result, which at the same time corrects the proof in [3], which went astray 
by an incorrect application of a time inversion argument, as explained in section 2. 



2 Expectation and variance of the area 

In this section, we prove the following theorem. 

Theorem 2 Let, for < a < b < oo, D{a, b) be the area of the region, bounded on the right and left 
by vertical lines through the x-coordinates of the points W{a) and W{b), and bounded from below 
and above by the line y = and the (left lower) boundary of the convex hull of V , respectively. 
Morever, let, for < a < b < oo, j3 = b/a and a = P — 1. Then 

(i) ED{a,b) = \\ogP, 

(ii) 

Var(D(a,6)) 

14, ^ 2 4 44 2{3 + a(3-4a)}tan-i(V^) 4. _w ^nn2 
= 27^°^^+3^ + 9^-45--^ ^ + 9(tan ^V^)) • 



Proof. We follow the line of the argument in Appendix A of [3]. First of all, we get, as in (4.3) of 
[3], replacing A{a,h) by D{a,b), 



and hence 



]im h-^ E {D {a, a + h)\W (a)} = ^aV{a f, 



/a pa 
sEV{sf ds = ^ J s-^EVil)"^ ds = l log a. 



Stationarity in the logarithmic scale then also gives ED{a, 6) = i log(6/a) and (i) follows. 

To prove (ii), we start with relation (4.4) in [3], which is repeated below (with A{l,a) replaced 
hjD{l,a)y. 

ED(l, of = ^ s'^EV{sfds + ^ E {D{1, s)V{sf} ds. (2) 

As noted on the top of page 46 of [3], the first integral is equal to (22/35) log a. As indicated on 
page 46 of [3], the other integral can be computed by a time reversal argument. But this is also 
exactly the point where things went astray. 

We have to compute the expectation of the product D{l,s)V{s)'^. D(l,s) is area of a con- 
catenation of trapezia, as for example shown in Figure 1. Note that Figure 1 corresponds to 
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w(i) = ([/(!), ni)) 




Figure 1: D{l,s) = area of shaded region 

three successive values of the jump process s D{1,s). For reasons of symmetry it is clear 
that ED{l,s)V{s)^ will be equal to ED'{l/s,l)U{l/s)'^, where D'{a,b) is the area of the region, 
bounded from above and below by horizontal lines through the y-coordinates of the points W{a) 
and W{b), and bounded on the left and right by the y-axis and the (left lower) boundary of the 
convex hull of V, respectively, sec Figure 2. The advantage of working with D'{l/s, in- 
stead of D{1, s)V{s)^ is that the somewhat elusive concatenation of trapezia, corresponding to the 
"developing area" (to use a phrase, coined in [3]), now lies in the future of the process with respect 
to U{l/s), whereas otherwise, in working directly with D(l, s)V{s)*, we would have to condition 
on the future value V{s) of the process. 
We have: 

limh'~^E{D'{a,a + h)\W{a) = 
ry 

= u (ux + ^au^) du 
Jo 

= + 

The situation is illustrated in Figure 3, where, during the time interval [a,a + h], a new point 
W{a + h) is caught, leading to the increment D'{a, a + h), which is a trapezium, bounded on the 
left by the y-axis. If W{a) = (x, y), the new point W{a + h) has coordinates (x + (a + h')u, y — u), 
where h' G [0, h] and u G (0,2/). 

Hence, using the time inversed process a i->- R[a) = {S{a),T{a)) = {V{l/a),U{l/a)), as on 
page 343 of [5], we get: 

E{D{1, a)\Wia)} = E{D'{l/a, l)|W^(l/a)} 

= f {lE{S{s)T{sf I R{l/a)} + lE{sT{sr | R{l/a)}} ds 

Jl/a 

s-mE{S{l/s)T{l/sf I R{l/a)} + lE{s-^T{l/sf \ R{l/a)}} ds 
= £ {l^^ {Vis)U{sf I Wia)} + \E {U{sf \ W{a)]] ds, 

and, instead of {5 / 12) EU{s)'^V{a)'^, wc have to compute: 

yEU{sfV{s)V{af + lEU{s)^V{a)\ (3) 

Note that the mistake in computing the second integral in (4.4) of [3] was caused by not taking 
into account that the relevant trapezium with upper right vertex W{l/s) is bounded on the left 
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W{l/s) = {U{l/s),V{l/s)) 




D'{l/s,l) 



W{l) = {U{l),V{l)) 



(0,0) 

Figure 2: £)'(l/s, 1) = area of shaded region 



D'{a,a + h) 




W{a) = {x,y) 

W{a + h) = {x + {a + h')u, y - u) 



x + ay x+ {a + h)y 
Figure 3: The increment D'{a, a + h) in the time interval [a, a + h], Note: h! G [0, h\. 



side by the y-axis and bounded below by a horizontal line above the x-axis, as is shown (for several 
of such trapezia) in Figure 2, instead of being bounded below by the rc-axis and bounded on the 
left by a vertical line to the right of the y-axis as in Figure 1. 

As in [3] , we can use the stationarity in the logarithmic scale to reduce the computation of (3) 
to the computation of 

\EU{lfV{l)V{sf + lEU{lYV{s)\ s > 1. (4) 
This expectation is computed in the the appendix. The result is: 

lEUilfvmis)^ + lEUiirVis)^ = 210 + M85 + 4(.-l).} _ 2(7 + .)^tan-M V^) 



15cr4(l + o-) 



where a = s — 1. We note in passing that the right side of (5) tends to 104/315, as cr I 0, so there 
is no singularity at cr = 0, as is also clear from the expression at the left side of (5). Let 

(Pis) = lEU{lfV{l)V{sf + lEU{lfV{s)\ s > 1. 

Then we can write, using (5), 

, 24 + 2s{105 + 4s(s-4)} 2(6 + s) tan'^ ( ^/J^) 

(P{S) = TZ—, — — : TTTTTT^ , S > 1. 



15s(s - 1)^ 



[S - 1)9/2 
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Thus the computation of the second integral in (2) boils down to the computation of 



We have: 



^ 4-19s + 2(s-l)3logs {s(14 + 5s) - 4}tan-i (Va) 104 



J u(f){u) du = 



^ / u<p(u) du = ^\ — h 



9o-3 9a7/2 945 

Thus 



5 rds r ^..^ 

jK I — / uMu) du 
Ji s Ji 



= i(loga)2- i24loga + | (tan-i (V^))' 

2{4 + a(4a - 11)} tan"! (Va^) 2 + 4a 44 
^ 9(a - 1)5/2 ^9(a-l)2 ~ 45 ' 

Using the notation a = a — 1, this can also be written 

— j u(t){u) du 

-in^^„\2 104,^^„ , 4 /.^„-w /-xx2 2{3 + Q(3-4a)}tan-^(V^) 6 + 4a 44 
-9(log«) - 945 log«+ 9 (tan (Vajj + ~9^ ~ 45 ' 

Hence, taking into account that the first term on the right side of (2) equals (22/35) log a, we get 
that the right side of (2) is equal to: 

1/, x2 , 14, , 4 -w ^\n2 2{3 + a(3-4a)}tan-i(Va) ,6 + 4a 44 
l(loga)^ + # loga + I (tan ^ (V^)) - ^ + " ' 

□ 



3 Appendix 

Proof of (5): By (4.6) in [3] we have, as a consequence of Lemma 2.4 on page 339 of [5], 
EU{l)^V{sf =lEUil)^ |F(l)4exp{-lc7F(l)2} - ^ exp {-1^^(1)2} 

+ ^{l-exp{-kni)n}}, 
where a = s — 1. Again using Lemma 2.4 in [5] we have: 

EU{l)^exp{-laV{lf}= [ x"^ exp {-^ay^ - l{x + yf} dxdy. 
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Let the integral Jj^2 a;^ exp { ^ay"^ — ^(x + y)^} dxdy be denoted by (f>{a). Then 

3(1 + a)2 tan-i (y^) 3 + 5a 



(^) 



r5/2 0-2 



This can be verified by using the recursive relations of Lemma 4.1 in combination with (4.8) in 
Note that 

0-4.0 o 

Furthermore, 

EUdfVdf exp } = = 3(l + '^)(S + <;)t'°--(v^) _ ]^ , 

and 

EU(lfV(lfe.,{-^,.V(lf}=M» = 3{3^ + 3.(10 + .)} tan- (y/?) _ ^21 + 11.) 

This yields: 



which coincides with (4.9) in [3]. 
For the other term we have: 



EU{lfV{l)V{s)'^ 

= lEuiif {F(i)^exp{-i.y(i)^} - yj^r^M-i^vm 

2F(l)(l-exp{-lani)'}) 



+ 



(72 



where a = s — 1. We have: 

EU(ifvmeM-hnm = I + ^ - ^"^tv"'''^' 

Define V(s) = ^C/(l)3y(l) exp {-ic7y(l)2}. Then 

£;t/(l)»F(l)3exp{-l.F(l)^} = -2/(.) = - 3(5 + 3-^) t^^"' (V^) _ 

and 

Eu(ifv(if..,{-i.vm = 4^"(.) = ^"^t^a+.V""' - 

Hence 

£;[/(i)V(i)y(s)^ 

_ 8 {210 + a(215 + 2(7(11 + a))} 24(14 + 5a) tan-^ (V^^) 16 

~ 25(t4(1 + (t) 5^975 ^ ^^4 0. 

Formula (5) now follows. 
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4 Concluding remarks 



There is a remarkable analogy between the behavior of the left-lower convex hull of the Poisson 
point process, discussed above, and the least concave majorant of (one-sided) Brownian motion 
without drift, as analyzed in [4]. In the same way there is an analogy between the behavior of the 

lower convex hull of the Poisson point process inside a parabola, as analyzed in [5] and [10], and 
the least concave majorant of Brownian motion with a parabolic drift, as analyzed in [6] and [8]. 
Why this is the case is still somewhat of a mystery and deserves (in my view) further investigation. 

Acknowledgements I want to thank Tomasz Schreiber for sending me the unpublished preprint 
[9] and Christian Buchta for pointing out some typographical errors in the original version of this 
manuscript (written in 2006). 
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